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1. Introduction 

Effective field theories are routinely used in particle physics and have proved an in- 
valuable tool for computing physical quantities in theories with disparate energy scales 
(for recent reviews see ||). Several years ago, Weinberg proposed that the machinery of 
effective field theory (EFT) could be applied fruitfully to nucleon-nucleon scattering and 
nuclear physics 0. Nucleon interactions might be profitably treated by EFT since they 
involve several different physical scales, such as the nucleon mass M, the pion and vector 
meson masses (m n , m p , m u , etc.). Furthermore, chiral symmetry in nucleon-pion inter- 
actions is necessarily expressed in the language of EFT, and the chiral expansion around 
m n = gives one a natural expansion parameter. Since Weinberg's original papers, much 
work has been done in the subject, with fair success in reproducing low energy features of 
nucleon-nucleon scattering from a chiral Lagrangian description of nucleon interactions || 

The goal of an EFT description of nuclear physics is not to improve upon semi- 
phenomenological models of the nucleon-nucleon interaction, such as the Paris @ , Bonn 
0] or Nijmegen [|J potentials . Instead it has been used to relate 3-body forces to 2-body 
forces H P] and to explain the observed hierarchy of isospin violation 0. One can also 



investigate the role of strangeness in hypernuclei or dense matter, along the lines of ||10|| . 
More generally, it allows one to better understand the physical origin of various features 
of the nucleon interaction (for recent progress in this direction see |TT[ - |T^] and references 
therein). One may also hope that the technique will allow semi-analytical approaches to 
solving 2- and many-body problems now only approached numerically. 

A fundamental difficulty in an EFT description of nuclear forces is that they are nec- 
essarily nonperturbative, so that an infinite series of Feynman diagrams must be summed. 
Which diagrams must be summed is well known, and the summing them is equivalent to 
solving a Schrodinger equation. However, an EFT yields graphs which require renormal- 
ization, giving rise to a Schrodinger potential which is too singular to solve conventionally. 
In Weinberg's work 0, only a contact interaction was summed, and the system was renor- 
malized; in @ , more complicated interactions are considered and a momentum cut-off 
is implemented, with bare couplings chosen to best fit phase shift data. In this paper we 
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focus on the 1 S (np) partial wave, and show how to compute the phase shift beyond lowest 
order in the EFT expansion, using dimensional regularization and the minimal subtraction 
(MS) renormalization scheme. 

Another problem with discussing systems with barely bound (or nearly bound) states 
in the language of EFT is that a new length scale emerges that is not directly associated 
with any physical threshold — the scattering length a. This makes the power counting in 
an EFT with large scattering length much less obvious than in one without. 1 5'o nucleon- 
nucleon scattering is particularly problematic from an effective field theory point of view, 
since the scattering length is very large: a ~ —24 fm ~ (8.5 MeV) -1 , a mass scale far 
lower than any hadron mass. In this paper we propose a specific ordering of the EFT 
expansion to avoid this problem. In the process, we are led to a modification of the power 
counting scheme proposed by Weinberg. 

We begin by briefly reviewing Weinberg's power counting scheme and the connection 
between Feynman diagrams and the Schrodinger equation. We then show how to sum 
the relevant graphs even when they are divergent, and we construct the low energy EFT 
for nucleons alone in the MS scheme. Finally we construct the EFT including one pion 
exchange in the MS scheme; in both cases we show at what scale the EFT fails. We 
conclude with thoughts about improving the approach, and its applicability to finite density 
calculations. 

2. Effective field theory, power counting and the Schrodinger equation 

2.1. Weinberg's power counting scheme 

The philosophy of EFT is that for scattering processes involving external momenta 
< Q, one need only consider a Lagrangian which explicitly includes light degrees of freedom 
for which m < Q. The effects of heavy virtual particles appear as an infinite number of 
nonrenormalizable operators suppressed by powers of the mass scale A relevant to the 
degrees of freedom excluded from the theory. EFT's can be predictive since amplitudes 
may be expanded in powers of Q/A, so that the effect of a nonrenormalizable operator on 
low energy physics is less important the higher the dimension of that operator. 
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The scale A can be determined by fitting low energy data to the predictions of the 
EFT to sufficient accuracy; the lower the scale A, the smaller the momentum range over 
which the EFT is predictive. An EFT will have to be modified as one approaches Q ~ A, 
and the degrees of freedom with mass A must then be explicitly included in the EFT. 
One then has a new EFT characterized by a scale A', which characterizes the next level 
of particles excluded from the theory. An EFT is only useful to the extent that there is a 
well defined hierarchy of mass scales; if there is such a hierarchy one can typically predict 
a large amount of data in terms of a few parameters. 

A necessary ingredient for an EFT is a power counting scheme that tells one what 
graphs to compute to any order in the momentum expansion. We reproduce here Wein- 
berg's analysis for NN scattering, couched however in the language of covariant rather 
than time ordered perturbation theory. The main complication arises from the fact that a 
nucleon propagator S(q) = i/(qo — q 2 /2M) scales like 1/Q if qo scales like m n or an exter- 
nal 3-momentum, while S(q) ~ M/Q 2 if qo scales like an external kinetic energy. Similarly, 
in loops f dqo can scale like Q or Q 2 /M, depending on which type of pole is picked up. 
To distinguish between these two scaling properties we begin by defining generalized "n- 
nucleon potentials" comprised of those parts of connected Feynman diagrams with In 
external nucleon lines that have no powers of M in their scaling 0. Such a diagram always 
has exactly n nucleon lines running through it, since there is no nucleon-antinucleon pair 
creation in the effective theory. includes (i) diagrams which are n-nucleon irreducible; 
(ii) parts of diagrams which are 1-nucleon irreduciblei. To compute the latter contribution 
to one identifies all combinations of two or more internal nucleon lines that can be 
simultaneously on-shell, and excludes their pole contributions when performing the J dqo 
loop integrations. An example of the 2-pion exchange contributions to is shown in 
Fig. 1. 

A special comment must be made about the 1-nucleon potentials, These dia- 

grams consist solely of the 1-nucleon irreducible graphs. They include both wave function 

1 With the exception of inverse powers of M from relativistic corrections. 

2 An n-nucleon irreducible diagram is one which does not fall apart when n nucleon lines are 
cut. 
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i<%. i. One /oop, 2-pion exchange Feynman graphs which contribute to the 
2-nucleon potential V^ 2 \ The first four are 2-nucleon irreducible; the last 
diagram is 2-nucleon reducible, and the poles from the slashed propagators 
are not included in the J dqo loop integration. The 1-loop graphs correspond- 
ing to tiNN vertex renormalization and pion wave function renormalization 
are not pictured here, but enter at the same order (as does nucleon wave 
function renormalization). 

renormalization (which begins at order Q 2 ) as well as relativistic corrections to the nu- 
cleon propagator, which start at order Q 4 /M 3 . The structure of the latter terms is fixed 
by relativistic invariance to reproduce the Taylor expansion of \/p 2 + M 2 . 

A general n-nucleon Feynman diagram in the EFT can be constructed by sewing 
together the nucleon legs of potentials with r <n; one treats the V^^s like vertices 
and the J dqo loop integrations pick up the poles of all the connecting nucleon lines!. 

The reason for this construction is that within the potentials, all nucleon prop- 
agators are off-shell and scale like 1/qo ~ I n contrast, when one picks up the pole 
contribution from one of the nucleon lines connecting the "vertices", other nucleon 
lines will be almost on-shell, and scale like 1/{Q 2 /M). 

Following Weinberg's arguments @, a contribution to the r-nucleon potential 
with I loops, I n nucleon propagators, I n pion propagators, and Vi vertices involving rij 

3 As pointed out by Weinberg, this set of diagrams is more naturally described in the language 
of time-ordered perturbation theory, but as there will be a mix of relativistic pion propagators 
and nonrelativistic nucleon propagators, no formalism is ideal, and we will keep to the language 
of covariant perturbation theory. 
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nucleon lines and d{ derivatives, scales like Q M , where 



\i = Al - I n - 2I n + J2 V i d i > ( 2 - la ) 
£ = I n + I K -^2Vi + l, (2.1b) 
In + r = l^Vim . (2.1c) 

In this power counting we take ~ Q and treat factors of the u and d quark masses 
at the vertices as order Q 2 . Combining these relations leads to the scaling law for the 
r-nucleon potential (r > 2): 

fi = 2 + 2*-r + £V i (d i + §n<-2) . (2.2) 

i 

Since chiral symmetry implies that the pion is derivatively coupled, it follows that (di + 
T^Tii — 2) > 0. That implies that for a 2-nucleon potential, \x > 0, and that fx = corresponds 
to tree diagrams. 

It is straight forward to find the scaling property for a general Feynman amplitude, 
by repeating the analysis that leads eq. fl2.2j) , treating the potentials as r-nucleon 
vertices with \i derivatives, \i given by eq. (|2.2|). However, while eq. ( |2.2|) was derived 
assuming that f dqo ~ Q and nucleon propagators scaled like ~ 1/Q, we now take them to 
scale like Q 2 /M and 1/(Q 2 /M) respectively. A general Feynman diagram is constructed 
by stringing together r-nucleon potentials V^ r \ 

For two nucleon scattering the situation is particularly simple, since the diagrams 
are all ladder diagrams, with n insertions of V^'s acting as ladder rungs. Each loop of 
the ladder introduces a loop integration (dq , od 3 q ~ Q 5 /M) and two nucleon propagators 
(~ M 2 /Q 4 ) to give a net factor of (QM) per loop. If we expand = J2^=o wnere 
V/2 ~ Q^, then a 2-nucleon diagram whose i th rung is the generalized potential V^f 
scales as 

L-l 

Q U (QM) L , v = > ( 2 ~ bod y scattering) (2.3) 

j=l 

where L is the number of loops (external to the V^'s). With insertions of V^- 1 ' along the 
nucleon propagators, which serve as relativistic corrections, there will be additional powers 



Fig. 2. The first two terms in the EFT expansion for the Feynman ampli- 
tude (T -matrix) for nucleon-nucleon scattering in the center of mass frame. 
The leading amplitude Aq scales like Q° and consists of the sum of ladder 

(2) 

diagrams with the leading (n = 0) 2-nucleon potential Vg at every rung; 
the subleading amplitude A\ scales like Q; it consists one insertion of 
(1-loop nucleon wavefunction renormalization) or one insertion of the sub- 

(2) 

leading (\x = \) 2-nucleon potential , dressed by all powers of the leading 
interaction Vq 2 ^ . 

of (Q 2 /M 2 ); likewise, an expansion of retardation effects in can be treated like the 
nonrelativistic expansion. 

Since fa > 0, the leading behaviour of the 2-nucleon amplitude is (QM) L . If one treats 
M ~ Q°, it follows that perturbation theory is adequate for describing the 2-nucleon system 
at low Q. In order to explain the nonperturbative effects one sees (e.g., the deuteron, or 
the large scattering length in the 1 5'o channel) one must conclude that M ~ 1/Q in a 
consistent power counting scheme. Thus the effective field theory calculation must be an 
expansion in v, given by eqs. Q2.2Q , (|0|). To leading order (v = 0) one must sum up all 

(2) 

ladder diagrams with insertions of Vq potentials with \i = 0. At subleading order one 

(21 (21 

includes one insertion of and all powers otVf 1 , etc. 

The program advocated by Weinberg is to solve the Schrodinger equation with the 
kernel expanded to a given order in \i. An alternative one might consider is to 
expand the Feynman amplitude A in powers of v; this is an equivalent procedure at v = 0, 
but is unsatisfactory for higher v as the expansion violates unitarity. We will argue in 
subsequent sections that for systems with a large scattering length (e.g. AW scattering) 
the best procedure is to expand |p|cot<5(p) = (47r/M)Re[l/*4] in powers of v, where p 
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is the momentum in the center of mass frame, and 5 is the phase shift. As discussed 
below, this expansion preserves unitarity and is expected to converge much faster than the 
expansion of the kernel V^ 2 \ particularly in systems with a large scattering length. 

A further modification we propose below is in the basic power counting scheme. In 
particular, we will argue that there are explicit powers of M in the coefficients of operators 
of the EFT, which decreases their order in the EFT expansion. We return to this point 
below, once we have presented the calculational techniques that lead to this result. 

2.2. Feynman diagrams and the Schrddinger equation 

Feynman diagrams are the usual tool for computing perturbative amplitudes requiring 
renormalization, while the Schrddinger equation is used to solve nonperturbative problems 
in potential scattering. As we will need to do both simultaneously, we briefly review here 
the connection between Feynman diagrams and the Schrddinger equation. 

Consider the integrals arising from the ladder loops in the diagrams of fig. 2: 

(2^)* V ( ^ ^' ^ (E/2 + q - q 2 /2M + it) {E/2 - q - q 2 /2M + ie) ^ ^' P ^ ' 

(2.4) 

In the above expression, M is the nucleon mass, E is the center of mass kinetic energy. 
(As we will focus entirely on the 2-nucleon problem for the rest of the paper, we will 
henceforth refer to the 2-nucleon potential simply as V). Following the rules of the 
previous section, the J dqo integral only picks up the pole contribution from the nucleon 
propagators at q = ±(E/2 — q 2 /2M). Since q ~ Q 2 /M one can consistently take q ~ 
in V (i.e., ignore retardation) to the order we will be working. In this approximation 

f -/(S y '"^) ( g - q VM + ^) y(q ' Pr (2 ' 5) 
The connection between the above expression and the Schrddinger equation is clarified by 

defining the free retarded Schrddinger Green's function operator for the 2-nucleon system 

G° E = 1 , H = — , (2.6) 

(E — H + ie) M 1 ' 

Matrix elements of G° E and the potential operator, V, between momentum eigenstates are 

given by 

= <Pmp'> = nP,P'). (2.7) 
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The sum of ladder diagrams can then be expressed as 

iA = -i ( P | (v + VG° E V + V(G° E V) 2 + ...) |p') 

= -i(p\V(l + G E V)\p') (2.8) 

= -z(p\(G E )- 1 G E (G E )- 1 \p'} , 
where G E is the full Green's function with potential V: 

Ge = ~ ^ ~ , H = Ho + V . (2.9) 
(E-H + te) 

We can define the state 

|X P ) = (1 + G E V) \p) = G^Gl)- 1 \p) (2.10) 

with p 2 /M = E, which is seen to be the exact scattering solution of interest: it satisfies 
the full Schrodinger equation since 

(H - E) | X p) = ~G E X | X p) = (Hq - E) |p) = , (2.11) 

and takes the appropriate asymptotic form for large r 

X P (r)^e l P- r + ^^e^, (2.12) 



since (r\G E V \r') oc 1/r for large r. The Feynman amplitude (|2.8| ) can be expressed in 
terms of x as 

iA = -i [ d 3 re- ip - r y(r)x P '(r) , (2.13) 



which is (— i) times the conventional expression for the T-matrix (see, for example, ||15| ). 

For s-wave scattering at center of mass momentum p, A can be conveniently related to 

the phase shift S(p) by the relation 

c/ \ ., , 4tt 1 
|p|cot(5(p) = i\p\ + — — 

;| M A (2.14) 

= — + lr p 2 + . . . , 

a z 

where a is the scattering length and r is the effective range. For 1 S'o (np) scattering, these 
parameters are measured to be (see |T 



a = -23.714 ±0.013 fm , r = 2.73 ± 0.03 fm . (2.15) 



The above discussion is complete when the potential V is less singular than 1/r 2 at 
the origin, in which case the terms in the series of eq. ( |2.8| ) are well defined. However, in 
effective field theories, the potential will in general have more singular behavior, such as 
1/r 2 , 1/r 3 , 5 3 (r), and worse. Such potentials do not allow a conventional solution to the 
Schrodinger equation, or equivalently, lead to divergent diagrams in the field theory. In 
the field theory it is well known how to deal with divergences — one merely regulates the 
integrals and then renormalizes the couplings of the theory, absorbing terms that diverge 
as the cutoff is removed into the definitions of the renormalized couplings. When this is 
done, there is no cutoff dependence in the theory. In this paper, we show how to sum up 
the leading diagrams using dimensional regularization and the MS subtraction scheme for 
the case of 1 5'o nucleon-nucleon scattering. This is equivalent to solving the dimensionally 
regulated Schrodinger equation. The advantages of our procedure are that dimensional 
regularization with the MS scheme preserves chiral symmetry and simplifies computations. 
Since the renormalization scale \x introduced by MS (or any mass independent scheme, 
such as MS) only enters in logarithms, EFT power counting arguments are particularly 
transparent, unlike when a momentum cutoff procedure is used. 

3. The effective theory with nucleons alone 

Although the power counting of the previous section assumed Q ~ m n and explicitly 
included pion propagation, the analysis also applies to a lower Q regime where the pion 
plays no role. We analyze this case first as it is analytically more accessible and quite 
instructive. 

At very low energy A A scattering we may consider an effective field theory consisting 
solely of nucleon fields; all other degrees of freedom, such as 7r's, A's, p and u mesons 
have been integrated out, and their effects are subsumed in the coupling constants of the 
effective theory. Note that this effective theory has nothing to do with chiral symmetry; 
in fact, it treats the pion as very heavy compared to momenta of interest, which is the 
opposite of the chiral limit. 

The EFT consists of all local nucleon interactions allowed by rotational invariance, 
isospin symmetry (which we assume to be exact in this paper) and parity. For 2-nucleon 
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scattering, the only interactions that are of relevance are the operators with four nucleon 
fields, as well as relativistic corrections to the nucleon propagator; we will be able to ignore 
the latter to the order we are working. In such a theory, the only diagrams that contribute 
to the 2-body potential V are tree diagrams. It follows from eq. (|2.2| ) that each 4-nucleon 
operator has a scaling dimension \ii = di, where d% is the number of derivatives in the acting 
at the vertex. Eq. ( |2.3j ) then tells us that the leading contribution to the amplitude has 
v = 0, and that Aq is given by the bubble sum of 4-nucleon operators with no derivatives. 
At v — 2, A2 is given by one insertion of a 2-derivative, 4-nucleon operator, dressed by the 
no-derivative operator, as in fig. 2, etc. 

The effective Lagrangian for this theory is given by 

V 2 

2M 2 bK ' (3.1) 

- \C 2 (N^V 2 N) (N^N) + h.c. + ... 

where a are the Pauli matrices acting on spin indices, and the ellipses refer to additional 

4-nucleon operators involving two or more derivatives, as well as relativistic corrections 

to the propagator. The coefficients Cs and Ct of dimension (mass) are the couplings 

introduced by Weinberg 0; C 2 is a coupling of dimension (mass) 4 . The values of Cs, 

Ct, C<i are renormalization scheme dependent. 

Nucleon scattering in the 1 5'o channel only depends on Cs and Ct in the linear 

combination C = (Cs — 3Ct), and so the leading contribution to the potential is 

Vb(p,p , ) = C. (3.2) 

Similarly, on can show that while there are a number of operators with four nucleon 
fields and two derivatives, only the linear combination proportional to C2 in eq. ( |3.1| ) 
contributes to V% in the 1 5'o channel. It will be convenient for later discussion of the 
momentum expansion to define 

C2 = ^ , (3.3) 

where A is a parameter with dimension of mass. With this definition the next to leading 
order contribution to the 2-nucleon potential is 

V 2 M) = C\^±^\ . (3.4) 
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3.1. The v = calculation 

The ladder graphs in fig. 2 for the leading part of the amplitude Aq can be summed 



trivially with the kernel Vq in eq. (|3.2|) , since the expression (|2.8| ) is a geometric series. 



Using dimensional regularization one finds 

-iC _ -iC 

° " 1-CG°(0,0) " 1 + zCM|p|/4tt ' (3 - 5) 

The quantity 

G° (0,0) = / ! £^<q|6g>|q'>, (3.6) 

sv ' ; J (2tt)™ 7 (27r) nXMI E IM 7 v ; 

is simply the coordinate space representation of the free Green's function Gg(r, r') = 

(r| 1/(E — Hq + ie) \r') evaluated at r = r' = 0, with reduced mass M/2. This corresponds 

to a divergent one loop graph, which in dimensional regularization is given by: 



G° E (0,0) 



d n q 1 



(27r) n (E -q 2 /M + ie) 
= - (An)~ n/2 M(—ME - ze) (n " 2) / 2 r(l - n/2) (3.7) 
My/— ME - ie __ -iM|p| 

n^3 4-7T 4-7T 

Even though minimal subtraction introduces a renormalization scale /i, one finds that 
in dimensional regularization (3^,(0,0) is finite as n — > 3 and so C is independent of /i, 
satisfying the trivial renormalization group (RG) equation 

4- (^1=0- (3-8) 



dfx \C j 

The value of C is determined by experiment via eq. ( |2.14j ) which fixes the threshold 
amplitude to be A = —Ana/M, where a is the scattering length. It follows from eqs. 
and (ET7D , using M = 940MeV, that 

47ra_ / 1 \ 2 
M V25 MeV / 

The expression for the scattering amplitude (|3.5|) may be rewritten as 



4tt/M 

Uo = z - -, ' . I I , 3.10 
—1/a — z|p| 
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which is recognized as the effective range theory expression for the amplitude given a 
scattering length a and effective range r = 0: 

Iplcottffp) = -- . (3.11) 
a 



This is reasonable, since the interaction (|3.1| ) is local. In fact, it is shown in the appendix 



how ( |3.9p may be derived by solving the Schrodinger equation with a potential V(r) = 
C5^(r). 

3.2. The fi = 2 calculation 

When using an effective Lagrangian, it is important to know at what momentum it 
fails. The momentum expansion in the effective Lagrangian ( |3.1| ) breaks down at the scale 
Q ~ A, where A is the scale set by V2 ( |3.4j) . To determine A, we must perform a second 
order calculation, at v = 2. How to do so is ambiguous: if one expands the Feynman 
amplitude to order v = 2 as in fig. 2, one destroys unitarity. In this section we will follow 
Weinberg's prescription, namely to expand V to second order (i.e., fi = 2), and then sum 
its effects on the amplitude to all order. Doing so includes the exact expressions for the 
order v = and v = 2 parts of the full amplitude, and keeps parts of the higher order terms 
(from multiple insertions of V2). In the following section we will consider an alternative 
calculation. 

To next to leading order, the 2-nucleon potential V is given by 

V(p, p') = V + V 2 = C (l + p2 ^ 2 p/2 ) . (3.12) 



It is possible to sum all of the ladder diagrams with the vertex (|3.12|) in the MS scheme; 
as shown in appendix B, one merely replaces C in eq. ( |3.5| ) by (7(1 + p 2 /A 2 ): 

lAv2 = l/[C(l + p 2 /A 2 )]+zM|p|/47T ' (3-13) 
where the subscript V2 denotes that we have followed Weinberg's prescription and expanded 
V (rather than A) to subleading order \x = 2. Since E = p 2 /M, the above expression for 
the amplitude can be expressed as a prediction for |p| cot5(p) by means of eq. ( |2.14j ): 

|PW> = -(g) C(1 + p V AV < 3 ' 14 > 
12 



We can fit our two free parameters C and A 2 to low energy scattering data by expanding 
|p| cot<5(p) to order p 2 and fitting to the measured scattering length and effective range 
ffnSD - The result is 

C= W = -{WMev) ' A^ =2roa = -l35MevJ ' (3 ' 15) 

Since p ~ A is the scale at which Vo — V2, we expect the effective theory with nucleons 
alone to work well at center of mass momenta |p| <C 35 MeV, but to fail completely for 
|p| > 35 MeV, corresponding to the lab kinetic energy Tj a b = 2.6 MeV. 

3.3. An alternative: expanding |p| cot<5(p) to order v = 2 

The result ( |3.15| ) is very discouraging from the EFT point of view. The original 
premise in §2 was that amplitudes could be expanded in powers of (Q/A) 1 ', where A was 
a mass scale typical of the particles not included explicitly in the theory. When pions 
are included, we would hope that A ~ m p \ in the lower energy EFT we are considering 
here, with the pion integrated out, one would expect A ~ m^. Instead, eq. ( |3.15| ) has 
A ~ 1/y/aro; ro can be considered a relatively short distance scale, but a ~ —1/(8 MeV) 
for the 1 5'o channel, which can hardly be called a typical QCD scale. In general, a blows 
up as a bound state (or nearly bound state) approaches threshold. Thus a small change 
in short distance physics can make the EFT fail at arbitrarily low momenta. 

The problem can be made more precise by examining the quantity |p| cot<5(p). Eqs. 
( jnp and ( FTI5D imply that 

1 1 °° 

|p| cot J(p) = = — ^(-iar oP 2 ) n , (3.16) 

-a + ia 2 r p 2 a ^ 

which has a radius of convergence at p 2 ~ l/(aro). However, it is known from general 
arguments that for a potential that falls of exponentially as e~ mr for large r, the true 
radius of convergence for |p| cot5(p) is given by p 2 ~ m 2 [TJ]. The quantity |p| cot<5(p) 
should have an expansion of the form 

1 00 

|p|cot<J(p) ~ -- + |r oP 2 ^(r 2 p 2 )" , (3.17) 

n=0 
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where the scales r n are typical of the range of the potential, tq ~ r n ~ 1/m. None of the 
r n 's are expected to diverge as \a\ — > oo. 

How are we to reconcile eq. ( |3.17|) with our discovery in eq. ( |3.16| ) that the scale 
of momentum variation in the EFT is set by the length scale y/afo ? The only possible 
answer is that the higher derivative operators in the EFT, although controlled by a scale 
that diverges as \a\ — > oo, are actually highly correlated, and the effects that diverge with 
a cancel. To see how this works in the present theory, consider a different expansion than 
the one performed above: instead of expanding V to order \i = 2 and solving for the 
amplitude, we will expand |p| cot<5(p) to order v = 2. In terms of a v expansion of the 
Feynman amplitude (as in fig. 2) 



A 



Ao + Ax + A 2 + 



(3.18) 



the expansion of |p| cot<5(p) is given by 
|p| cotJ(p) 



^ + MA 



z p + 



4tt 1 

MA 
4tt 1 

MA 



AA ( A\- A A 2 
Ao) + V Al 



+ . . 



(3.19) 



Note that to compute |p| cot<5(p) to order z/ , one needs to compute A u only for v < Vq, 
which involves perturbation theory in all but the v = potential @. In the present EFT, 
we have found 

Ait/M 



Ao 



C 



l + iCM\p\/4n 
(from eqs. flOD , (l3T5|) ), Ai = 0, and 



A 2 



C 



T/a — i|p| 



(3.20) 



l + iCM\p\/4n J \CA 2 J 



(3.21) 



obtained by expanding eq. (|3.13j ) to first order in 1/ A 2 and substituting the values ( |3.15|) . 
Substituting the above expressions into eq. ( |3.19| ), we find that the v = 2 expansion of 



4 In a theory with just nucleons, all A u vanish for odd v\ this is no longer true when pions are 
included. 
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Fig. 3. 1 Sq np phase shifts in degrees plotted versus center of mass momen- 
tum. The dots are the 1 Sq phase shift data from the Nijmegen partial wave 
analysis the dashed, dash-dot and solid lines are EFT calculations in 



a theory without pions. The dashed line is the v = result from eq. ( \3. 
the dash-dot line is the EFT result when the potential is expanded to order 
v = 2, eq. (\3.16j ; the solid line (which lies along the dots) is the EFT result 
when |p| cotJ(p) is expanded to order v = 2, eq. ( \3.2Q ). 

|p| cot5(p) exactly reproduces effective range theory, 

|p|cot(5(p) = — + ±r p 2 (^ = 2). 



(3.22) 



In fact, with no long range interactions, effective range theory had to be equivalent 
to the v = 2 EFT expansion of |p| cot5(p), by dimensional analysis. However, when pions 
are included the EFT expansion of |p| cot5(p) is not equivalent to effective range theory, 
since each order in the v expansion generates a complicated dependence for |p| cot5(p) on 
p and m-n- of order Q v . In fact, as we show in the next section where we include pions, the 
lowest order, v = calculation (with the scattering length a as experimental input) allows 
us to predict a nonzero value for the 1 5'o effective range r . 

3.4- Comparison with data 

In fig. 3 we show a plot of the 1 S'o iViV phase shift produced by the Nijmegen partial 
wave analysis ||17|| , compared with the three EFT analyses we have performed without 
pions: the v = calculation (|3.11| ); the /i = 2 expansion of the kernel V ( |3.16| ); and the 
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v = 2 expansion of |p| cot<5(p) (|3.22| ). The first two calculations yield results that agree 
well with data for |p| < 20 MeV, but differ wildly above that scale. This is what one would 
expect from the size of C and A 2 in eqs. Q3.9|) , (|3.15|) . In contrast, the v = 2 expansion of 
|p| cot<5(p) yields a result indistinguishable from the Nijmegen analysis beyond |p| = m^. 
This demonstrates the nearly complete cancellation between operators of different index 
|U discussed above, at momenta greater than A. 

We do not wish to belabor this phenomenological success with this simplistic model 
- we have only shown that low energy EFT can reproduce effective range theory, and 
we have only considered a single partial wave. However, we have demonstrated that ex- 
pansion of |p| cot<5(p) extends the range of validity of the EFT beyond the scale set by 
the derivative expansioni. There is also an important practical reason for preferring to 
expand |p| cot<5(p): that is that the effects of \x > interactions need only be computed 
in perturbation theory, following the expansion (|3.19 ). This in general leads to a great 



simplification of the calculation. Furthermore, it provides a way to implement a consistent 
renormalization procedure, as we will see in the next section, where we introduce pions. 

3.5. Rethinking the EFT expansion 

We showed above that the effective Lagrangian is a momentum expansion in powers 
of p 2 aro, where r$ is a typical hadronic scale (ie.~ 1 fermi), while a is the scattering 
length. This led us to the conclusion that the object with a sensible momentum expansion 
is not the potential, but rather Re [1/-4], the real part of the inverse Feynman amplitude. 
Another conclusion we can draw is that the power counting scheme presented in §2 needs 
modification — eq. ( |2.2|) in particular. The leading, 4-fermion operator has a coefficient 
4-Ka/M, and is treated as a \i = contribution to the potential. Since the combination MQ 
is assumed to be of degree v = 0, it follows that powers of aQ appearing in the amplitude 
are also of degree v = 0. As the effective Lagrangian is an expansion in ~ p 2 aro, a 2d 
derivative vertex does not contribute \i = 2d to a graph, as assumed in eq. (|2.2| ), but 
rather n = d as it scales like (aQ 2 ) d . This does not matter in the present theory; the 



5 This suggests that the EFT might be profitably formulated with a light degree of freedom in 
the s-channel. 
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u u = 2 calculation" needs only to be renamed the u u = 1 calculation". However as we 
show below, there are striking implications when pions are included in the EFT. 



4. The effective theory with nucleons and pions 

In order to extend the energy range over which the effective theory is useful, it is 
necessary to include more light degrees of freedom. The obvious candidate to add to the 
theory is the pion. In the previous section, analyzing np scattering in an effective theory 
without pions, we found that the derivative expansion in the EFT broke down at a scale 
A ~ 35 MeV. A sign that we are improving the utility of the EFT by including pions will 
be whether or not the scale set by the contact interactions becomes significantly higher. 
As we will show, that is the case. 

Chiral symmetry mandates that pions couple to nucleons derivatively, or proportional 
to powers of the quark masses. In the power counting arguments of §2, we assumed 
Q 2 ~ m 2 oc m qi where m q are the u and d quark masses. To determine which operators 
to include at a given order in the EFT expansion, it is necessary to look to eqs. ( |2.2j ) 
and (pT3|). To compute the 2-nucleon potential at order \i = we include all tree level 
interactions for which ^Vi(di + \ni — 2) = 0. That includes the 4-nucleon interaction 
without derivatives, as well as 1-pion exchange with the 1-derivative axial vector coupling 
at each vertex. In this section we perform the leading [y = 0) and subleading {y = 1) 
calculations exactly. 

4-1. The v = amplitude 

The \x = 0, 2-nucleon potential Vo for NN scattering in the 1 5'o channel is given to 
leading order by one pion exchange, plus a contact term: 

v( P , P ') = c - (fi) (q -7r; 2) ^ i,T2> ■ 

with q = (p — p'). The coupling g& = 1-25 is the axial coupling constant, = 140MeV 
is the pion mass, and f n is the pion decay constant normalized to be 

U = 132 MeV , (4.2) 
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Fig. 4- Ladder diagrams for the v = contribution to the Feynman am- 
plitude Aq are formally resummed by expressing the kernel Vq as a sum of 
a contact interaction proportional to C and a nonlocal interaction V n as in 
eq. The shaded blobs consist of the ladder sum of V n interactions 

(dashed lines), while the black vertices correspond to a factor of C . 

compared to other common normalizations f n = \/2(93 MeV) = (186 MeV)/\/2- As in 



the previous section, C is a free parameter which will be computed in MS subject to the 
condition that we correctly reproduce the measured threshold scattering amplitude (i.e., 
the scattering length a). Since we are exclusively interested in the 1 S'o channel (I = 1) we 
can express Vq as 

Vb(p,p') = C + V 7t (p,p') , (4.3) 

where 

Note that while is the conventional one-pion exchange (OPE) potential, our calculation 
will differ significantly from OPE due to the C contact interaction. The contact term 
includes not only the 5 3 (r) contribution from one pion exchange, but also the leading 
contribution in the derivative expansion of all shorter distance effects, such as 2-pion 
exchange, intermediate A's, u exchange, etc. 

It is not possible to compute the ladder sum with the above kernel analytically, but 
we are able to express it in terms of several quantities that can be computed numerically 
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Fig. 5. Subdiagrams defining the quantities A n , X P (0) and Ge(0,0) used 
in eq. ( \4-W - *A n and x P (0) are finite, as are all but the first two diagrams 
in the expansion of Ge(0, 0) . The shaded blob is defined diagrammatically 
in fig. 4- Dashed lines are insertions ofV n , eq. ( \4-4i ) 

with ease. Most importantly, we are able to renormalize the nonperturbative amplitude 
analytically. To achieve this, the ladder diagrams are formally summed as in fig. 4 to yield 
the Feynman amplitude! 

iAo = iA n -i C[x g ( ° )]2 , (4.5) 
1-CG E (0,0) 

where the quantities A n , X P (0) and Ge(0, 0) are the sub-diagrams pictured in fig. 5. 

The quantity A n is just the amplitude one finds in the pure Yukawa theory with 

potential V n , i.e, the usual OPE result: 

iA !r = (p\%(l + G E %)\p') , 

1 (4-6) 



G 



E 



E-Ho-V^ + ie 
while 

d 3 q 



Xp(0) = J ^3 (ql C 1 + GeVx) |p) (4.7) 
is the OPE wave function at the origin. Both A n and x p (0) can be computed numerically 
by solving the Schrodinger equation with the Yukawa potential V n . This is discussed in 
appendix A, where the solutions are plotted (figs. 7,8). 

The quantity G\e(0, 0) is the coordinate-space propagator from the origin to the origin 



in the presence of V n ; it is divergent but can be defined in MS 



6 Here we give a Feynman diagram approach, while in appendix A we show how the Schrodinger 
equation corresponding to the kernel ( |4.1| ) can be solved directly. 
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Divergences occur in only the first two graphs in the perturbative expansion for Ge(0, 0) 



shown in fig 5. These two graphs can be computed in MS, while the remaining graphs 
can be summed by numerically computing the propagator Ge(t, 0); see appendix A for 
details. As a result of renormalization, Ge(0, 0) is replaced by the finite Gj^ s (0, 0), while 
the bare C is replaced by the renormalized Cj^-gip) which is to be fit to experiment. Note 
that both quantities now depend on a renormalization scale \x\ however the amplitude A§ 
is \i independent. We can compute the renormalization group equation for Cj^-g(fx), which 
is given by 

d 1 a n M< 



— = - A (MS) . (4.9) 

^« ^ 

This result is derived in appendix A (see eq. (A. 28)). Throughout this paper we will be 
quoting values for coupling constants renormalized at the scale \i = m^; the reason for 
this is that loop diagrams omitted at a given level of the v expansion bring in factors 
of lnm 2 /^ 2 , and so choosing [i ~ is expected to optimize the perturbation expan- 
sion for |p| < m n . Note that for Cjjg(p) negative, \1 / Cjjg(fi)\ increases with increasing 
renormalization scale \i. 

After solving for A n , x P and G^ s numerically, one can compute the amplitude 
and fit Cjjg(p) in order to obtain the correct scattering length. We find 



1 \ 2 



79 MeV 



(v = 0), (4.10) 



(Cj^c = Cjj-g — g 2 A /2f%) which shows a substantial improvement over the value C = 
— (1/25 MeV) 2 obtained in the pure nucleon effective theory at v = 0, eq. (|3.9j ). Fur- 
thermore, once C^-j-g is fixed to give the correct scattering length, one can compute the 
effective range, and we find 

r = 1.3 fm [y = 0) , (4.11) 

which shows that this simple effective theory with a contact term and one-pion exchange 
can account for about half of the measured effective range, tq = 2.7 fm. In fig. 7 we plot 
the phase shift determined from the amplitude Aq in eq. ( f4.5|) as a function of the center 
of mass momentum |p|. 
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(a) (b) 

Fig. 6. Examples of graphs with logarithmic divergences: (a) one propor- 
tional to M 2 m 2 ; (b) another with a term proportional to M 2 p 4 . The solid 
lines are nucleons, while the dashed lines are pions. 



4-2. A new power counting scheme 

In section §3.5 we pointed at that the parametrically large coefficients of operators 
in the effective Lagrangian forced us to revise our power counting scheme; however the 
revised counting had little practical effect in the theory without pions. However, with 
pions the story is different. In particular, equation ( |4.9| ) above mandates that we modify 
our EFT expansion scheme significantly. The logarithmically divergent graph that leads 
to the scaling of C is shown in fig. 6(a); in 4 — 2e dimensions it gives rise to the pole 

1 a n M 2 1 g\m\M 2 



16tt ~ e 128?r 2 /2 



(4.12) 



By our power counting, it is perfectly consistent to find such a divergence at v = 0, since 
m 2 M 2 has degree v = 0. However, note that the counterterm required to absorb this 
divergence is of the form (JVt Ai q N)(N^ N), where A4 q is the quark mass matrix. In the 
power counting scheme described in §2, such an operator, which is higher order in a chiral 
expansion, was considered to be a (jl = 2 operator, since it was assumed that the coefficient 
of the operator was set by a "typical" QCD scale. Instead, we see from eq. ( |4.12| ) that 



the coefficient of this operator has an explicit factor of M 2 , and so is actually of degree 
^ = 00. 



7 The connection between the coefficient of a 1/e pole, and the natural size of the coefficients 
of operators in the effective Lagrangian is "naive dimensional anaiysis". The idea is simpiy that, 
because of the RG flow such as in eq. (4.9), even if the coefficient C(^t) of an operator is small 



for some reason at a scale fio, at a scale u = uo X 0(1), C(/i) will have flowed to the magnitude 
of the coefficient of the 1/e pole, which is therefore considered its "natural" size. 
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Evidently the EFT expansion is not equivalent to a chiral expansion in m^, since a 
calculation at finite order in the EFT expansion can be modified by arbitrary powers of 
M 2 m 2 . This may sound like a disaster for the EFT expansion, but it actually is not. At 
v = 0, for example, we still have only one contact interaction C, and whether or not this 
contains contributions from all orders in has no effect on the predictive power of the 
calculation. 

What would be a disaster is if there were counterterms needed proportional to powers 
of M 2 p 2 . If that were the case, the theory would not be predictive, as an entire form factor 
would be needed to describe scattering at lowest order in the EFT expansion. However 
such terms do not arise. Consider the diagram fig. 6(b), which contributes to the v = 
2 calculation of the amplitude. The graph is proportional to C 2 M 2 /f 4 ~ l/(mass) 6 . 
Since the graph must be proportional to l/(mass) 2 , there can be a logarithmic divergence 
(and hence a 1/e pole that requires a counterterm) proportional to C 2 M 2 p 4 // 4 . This is 
consistent with the graph being v = 2; it is also consistent with the power counting of the 
theory without pions, where the p 4 interaction was multiplied by (aro) 2 . However, this 
result is not consistent with the naive power counting in §2 that assigned degree fx = 4 to 
a four nucleon contact interaction proportional to p 4 . 

It is possible to make general power counting arguments that the four nucleon contact 
terms in the effective Lagrangian involve an expansion in M 2 m 2 and Mp 2 . Equation ( |2.2j ) 
is therefore modified so that the degree fx of a vertex gets a contribution Afi = 1 rather 
than Afi = 2 from each factor of nucleon momentum squared, p 2 . Furthermore, powers 
of M. q in four nucleon contact interactions are assumed to be accompanied by M 2 and do 
not increase fx. 

4-3. The full v = 1 amplitude with one-pion exchange 

From the above discussion we see that the p 2 contact interaction — given that its 
coefficient scales ~ M — is the only operator in the effective Lagrangian of degree \x = 1. 
Two pion exchange — considered of equal order in the incorrect power counting of §2 - 
contributes at order fx = 2, as does the p 4 contact interaction. A full v = 1 treatment of 
the 1 Sq scattering amplitude is therefore simple to obtain, while extending the analysis to 
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v = 2 is quite ambitious. 

As in §3, we first consider the Weinberg expansion, summing up the \i 



V = Vq + Vi to all orders, where Vq is given in eq. (|4.3| ) and 



Vi = C 



P 2 + 

2A 2 



1 potential 



(4.13) 



The resultant amplitude is (see appendix B for details) 



%A 



[Xp(0) 



C (1 - a,m T M/A 2 + p 2 /A 2 ) 



G E (0,0) 



(4.14) 



where A w , X P (0) and (5^(0,0) are defined as in the previous section, fig. 5. Expanding 
the denominator of eq. (4.7) in powers of 1/A 2 gives 



1_ 1 / a^m^M p' 
C CV A 2 A 2 



G B (0,0) . 



(4.15) 



The term Ge is divergent and if it is defined using dimensional regularization it has an 



energy independent 1/e singularity. Using the MS subtraction scheme, this divergence 
is absorbed into a renormalization of C changing it to Cj^g{^). Then in the second 
term proportional to 1/A 2 we must introduce a renormalized A-^-^) defined by CA 2 = 
Cjjg((i)Ajj-g((i) 2 . However, with C and A renormalized in this way there is no freedom 
to express the higher order terms represented by the ellipses in eq. ( [4 . 1 5| ) in terms of 
renormalized parameters. This problem arises because we have not included operators 
with more than two derivatives. They are needed as counter terms to render multiple 
insertions of the two derivative operator in eq. ( |3.1|) finite. This is equivalent to saying 
that no redefinition of couplings in eq. ( [4.14j ) can absorb the energy independent 1/e pole 
in G E (0, 0). 

A procedure which we can follow, consistent to order v = 1 in an expansion of the 
amplitude, is to include all the higher derivative operators, absorb the 1/e, and then 
arbitrarily set the renormalized coefficients of the higher order terms to zero. This ad hoc 
procedure results in the analog of eq. (|3.13|) , with one pion exchange effects included: 

[Xp(0)] 2 



Cms ( 1 - a w m w M/A2_ + P 2 /Af_) 



MS 
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Gf s (0,0) 



(4.16) 
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Fig. 7. 1 Sq np phase shifts in degrees plotted versus center of mass mo- 
mentum. The dots are the 1 Sq phase shift data from the Nijmegen partial 
wave analysis j\ff[J; the dashed, dash-dot and solid lines are EFT calcula- 
tions in a theory with one pion exchange. The dashed line is the v = 
result from eq. ( U ■ o\ ) ; the dash-dot line is the EFT result when the potential 
is expanded to order v = 1, eq. (^..14 ); the solid line (which lies along the 
dots) is the EFT result when |p| cotJ(p) is expanded to order v = 1, eqs. 
( U-1% ) — ( U-2(\) . Note that the momentum range of the plot extends to 2m 7I -, 
twice the range of fig. 3. 

A (numerical) fit to the measured scattering length and effective range with the amplitude 
in eq. ( f4.16| ) gives 



C Ws^) 



1 



1 



A— (m) 

MS vr ' 



1 



(43 MeV) 2 



(4.17) 



(125 MeV) 2 ' 

Comparing the above values with the analogue without pions, eq. (|3.15|) , we see that 
the inclusion of pions has greatly reduced the value of C, while not significantly altering the 
scale A of the derivative expansion. Therefore, the range of utility of the EFT apparently 
remains disappointingly small. In particular, the scale A is far below the Fermi momentum 
in nuclear matter [pp ~ 280 MeV) so that the theory would appear to be of little utility 
in understanding nuclear physics. Nevertheless, following the discussion in §3.3, we expect 
an expansion of |p| cot5(p) to order v = 1 will work much better; this is indeed the case, 
as demonstrated in fig. 7. 

To compute the v = 1 expansion of |p|cot5(p) we need to know Aq and A\. The 
amplitude Aq was computed in the previous section in eq. ([4.5| ), with the substitution of 
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the MS values for C and Ge- 

A = A 



[Xp(0) 



(4.18) 



The second order amplitude A± corresponds to the sum of graphs indicated in fig. 2; their 
sum is given by our previous calculation ( [4 . 1 6| ) expanded to first order in 1/A 2 : 

[Xp(0)] 2 



At 



^ms 



(a^M/ A 2 ^ - P 2 M^) ■ (4-19) 



The v = 1 expansion of |p| cot5(p) is given by eq. ( |3.19| ): 



|p| cot 5(p) = i\p\ + 



4tt 1 
MA 



Ax 
A 



(4.20) 



This procedure is well defined from the point of view of renormalization: Note that aside 
from the explicit factor of p 2 in eq. ( |4 . 1 9|) , there is also complicated momentum dependence 
in An, X P (0) and G^f s (OPE Feynman amplitude, the OPE wave function at the origin, 
the renormalized OPE Green function at the origin respectively). Thus the terms in the 
expansion ( |4.2U[ ) do not correspond to the two parameters in effective range theory; indeed, 
we saw in the previous section that the v = contribution already accounts for half of the 
effective range. 

By fitting the two free parameters Cjjg and so that the expression ( |4.20|) cor- 
rectly reproduces the 1 5'o effective range and scattering length, we arrive at the prediction 
for the phase shift plotted as a solid line in fig. 7. The values one finds for the parameters 
are now: 



fi=m T , 



(121 MeV) 2 



(4.21) 



(100 MeV) 2 ' 

which indicates a very significant improvement over those found by first expanding the 
potential to order /i = 1 and then summing to all orders, eq. (|4.17| ). In particular, the 
momentum expansion scale A is now much larger. 

Even with the larger scale |A| = 121 MeV, one would not expect the momentum 
expansion to converge fast enough to be of use in nuclear matter, where pp ~ 280 MeV. 
However, as argued in the previous section and evidenced by fig. 7, the expansion for 
|p| cot5(p) has a much larger radius of convergence than the derivative expansion in the 
Lagrangian. 
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5. Conclusions 

We have shown how to perform a nonperturbative calculation of NN scattering in the 
1 5'o channel in an effective field theory expansion. A key feature of the procedure was the 
application of dimensional regularization (usually viewed as a perturbative regulator) and 
the MS renormalization scheme procedure to the nonperturbative problem. Our results 
for the phase shift depend only on physical observables and not on any momentum cutoff, 
even though the bare NN interactions in an EFT are singular. 

At leading order in the EFT expansion {y = 0), which includes one pion exchange and 
a contact interaction, we find a prediction for the effective range tq = 1.3 fm, given the 
measured scattering length; this is about half the measured value. The fit to the measured 
phase shift is poor above |p| ~ 25 MeV. In order to better understand the range of validity 
of the EFT approach, we investigated the phase shift including effects at subleading order 
in the EFT expansion. At this order there is an ambiguity about what quantity should 
be expanded; the ambiguity corresponds to which higher order terms are kept in the EFT 
expansion to maintain unitarity. Following the method of || , one can expand the potential 
to subleading order in the EFT expansion, and include its effects to all orders. Doing this, 
we find the phase shift that results disagrees with data above |p| ~ 45 MeV, which is 
what one expects from the size of coefficients one finds for the derivative expansion of the 
effective Lagrangian. 

An alternative method we explore is to expand the quantity |p| cot5(p) to subleading 
order. We explain why this expansion should be expected to have a greater radius of 
convergence than the derivative expansion would lead one to expect, at least at low orders 
in the EFT expansion. This is supported by calculation, which suggests that the 1 Sq phase 
shifts at subleading order agree well with data at up to ~ 280 MeV. A strong correlation 
is implied between coefficients in the derivative expansion of the Lagrangian that has to 
do with an s-channel pole, but which remains to be quantitatively understood. 

By investigating the EFT both with and without one pion exchange to order v = 1, 
we see that including the pion increases the inverse mass scales that appear in the EFT 
expansion, thereby improving improving its utility at high momentum. Including two 
pion exchange, four derivative interactions and possibly the effects of the A will increase 
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these scales even further. Since higher partial waves are less sensitive to short distance 
physics (and are in fact well approximated by one pion exchange, which appears at lowest 
order in the EFT expansion), we are optimistic that the techniques presented here will be 
successful at reproducing all of the spin singlet partial wave phase shifts up to center of 
mass momenta comparable to the Fermi momentum in nuclear matter. This investigation 
is in progress. 

As a consequence of our analysis we are forced to revise the power counting scheme 
of Weinberg || in order to account for the powers of the nucleon mass M that appear in 
operator coefficients of the effective Lagrangian. We find that the EFT expansion is not 
equivalent to the chiral expansion, as the EFT expansion requires summing contributions 
proportional to all powers of M 2 m 2 . 

Application of these techniques to the spin triplet channel is not straightforward, 
however, since the interactions in this channel are singular but not separable {e.g., a 1/r 3 
singularity from one pion exchange). Our hope is that this problem can be surmounted, in 
which case the techniques we developed here should prove of use in a variety of interesting 
problems. The EFT approach could be applied to nuclear matter, with the goal of un- 
derstanding its binding energy and compressibility in terms of a few parameters extracted 
from low energy scattering experiments. One could investigate the implications of 577(4) 
symmetry in N and A interactions, recently shown to be a consequence of the large- N c 
expansion of QCD |18|] [^9|. In particular, £77(4) symmetry greatly reduces the number of 



four-fermion operators one needs to consider when the A is included [19 . 



Since 577(3) flavor symmetry and its breaking can be easily incorporated in the EFT 
formalism, it may prove a useful tool for exploring systems with nonzero strangeness, 
extending the discussion of ref. |10| to a nonperturbative analysis. Finally, of great interest 
is the possibility that the EFT analysis may prove to be a useful tool in understanding 
systems at densities above nuclear density, with an eye toward a systematic inclusion of 
nuclear forces in the presently incomplete analyses of pion condensation [2(J and kaon 
condensation ]21[]-[^3[ . 
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Appendix A. Solving the Schrodinger equation 



In the §3 we showed how to sum ladder diagrams involving a 4-nucleon contact interac- 
tion. This is formally equivalent to solving the Schrodinger equation with a 5 3 (r) potential; 
however neither approach makes sense without renormalization. While the techniques of 
renormalization are familiar in the context of field theory, here we show how to obtain 



the same results via the Schrodinger equation (See also ref. f24 ]. This approach is quite 
convenient for practical computations. 
The equation we want to solve is 







where 



-V 2 /M + K-(r) + C5 3 (r) - E ^(r) 
H - E + CS 3 (r)] ifj(r) , 

VJr) = -o,,- . 



(A.l) 



(A.2) 



Away from r = we can find two independent s-wave solutions to (H — E)ip = We 
denote the regular s-wave solution by Je{t) and the irregular s-wave solution by JCg(r) 
They are normalized to have the following behaviour near r = 0: 

(A.3) 



Je{t) — > 1 



These functions have several features: 



4n 



In Ar + O(rlnr) . 
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(i) K,^ is a Green's function satisfying 

(H - E)K% = 5 3 (r) ; 

(ii) The arbitrary scale A in /C^ corresponds to the choice of boundary conditions on the 
Green's function (i.e, the arbitrariness in redefining K,e{t) by an amount proportional 
to J E (r))- 

(iii) For both functions, the dependence on the energy E vanishes as r — ► 0; 

(iv) asymptotically, these functions become: 

Je(t) — ► (ye ipr /pr + c.c.) , 

(A.4) 

4(r) — > (ze ipr /pr + c.c.) , 

r — >oo 

where E = p 2 /M and y and z are complex constants that must be determined nu- 
merically. However, y and z are related: since f [JC E (H — E)Je\ d 3 r = 0, it follows 
upon integration by parts that 



The Schrodinger equation ( |A.1|) can now be rewritten as 



(H - E)iI){t) = -CV(0)5 3 (r) , (A.6) 
and is formally solved in the s-wave channel by 

^(r) = aJ E (r) + blC E (r) — ► (ay + bz) + c.c. (A.7) 

r— >oo pr 

provided that b = — (7-0(0), or 

C 

b = -a ^- . (A.8) 

1 + C7C*(0) 

(We say "formally" since /Cg(0) is divergent; we will address the issue of renormalization 
below). Note that the ratio a/b is real. Comparing eq. ( |A.7| ) with the desired (s-wave) 
asymptotic boundary condition 

D ipr p — ipr ' 



2 

r — >oo 



if,(r) — > -| e M , (A.9) 
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50 100 150 200 250 300 

p cm (MeV) 

Fig. 8. The OPE phase shift 5 n (in degrees) as a function of |p| in MeV. 

it follows that the phase shift is given by 

e 2lS = -(ay + bz)/(ay + bz)* 



_y_ 

y* 



1 



zy — yz 



y* J a/b + z*/y' 1 



{ ipM 
\~8tT 



(A.10) 



XJ \ Syr J -l/C- }C E (0) + z*/y* 
where we have made use of eqs. (|A.5|) , ( |A.8|) , and have defined 5^ to be the "OPE" s-wave 
phase shift arising from the one pion exchange Yukawa interaction V n , and no contact term 
(exp(2^) = -y/y*) (see fig. 8) i. 

It is now just a few steps to relate the above expression to eq. (f4.5|), the analogous 
formula derived diagrammatically. First note that the canonically normalized scattering 
solution in the pure Yukawa theory is given by Xp( r ) = — iJE{f) / {2y*), so that 

Xp (0) = -i/(V) • (A.ll) 

The function [x P (0)] is plotted in fig. 9 as a function of the centre-of-mass momentum. 

Next note that the retarded Green's function (satisfying the asymptotic boundary 
condition that there is no incoming wave) is given by 

G E (r, 0) = (-/Ci(r) + (z*/y*)J E (r)) ; (A.12) 
8 Our numerical calculations were performed with = 140 MeV and M = 940 MeV. 
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50 100 150 200 250 300 

p cm (MeV) 

Fig. 9. The OPE wave function squared at the origin, [x P (0)] 2 , plotted 
versus momentum in MeV. The solid and dashed lines correspond to the 
real and imaginary parts respectively. 



finally, the relation between the Feynman amplitude iA and the phase shift is 



iA 



.47r(e 



2 id' 



2ipM 



(A.13) 



It follows that expression ( A.10 ) is equivalent to eq. (4.5) 



c[x P (o) 



1 - CG E (0,0) 



(A.14) 



When the effects of pions are not included, as in §3, one recovers the amplitude ( |3.10|) , 
since a n — ► implies A n — > 0, C — > C, X P (0) — > 1, and Ge — ^ G E in the above expression. 

So far the discussion has been in terms of the quantity, K, E (0), which was seen in eq. 
( |A.3| ) to have both linear and logarithmic divergences as r — > 0. This can be remedied by 
renormalizing C, for example by defining 



1 



i + /c A (o) = i + /c£(o) 



(A.15) 



c R (\) 

(where we used the fact that lim r ^ [)C E (r) -K% (r)] = 0, from eq. (gjg)). Then eq. (|OC| ) 
for the phase shift can be written in terms of renormalized quantities as 

1\ 2 fipM\ 1 
V*J V^J -l/C R (\)+z*/y* 
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2iS 



(A.16) 



This choice of renormalization scheme is convenient for computations: First one solves 
for the solutions to the Schrddinger equation with the OPE potential V n (r), subject to 
the boundary condition (|A.3|) ; from that one computes the asymptotic behaviour and the 
coefficients y and z in eq. (|A.4[ ); then one computes Cr(X) from eq. ( |A.16| ) in terms of the 
measured 1 So scattering length a = —23.7 fm. Note that this expression is very different 
than the pure OPE calculation. 



The renormalization prescription ( |A.15|) is different than MS, but it is straightforward 
to relate the two. Using 



one finds 



1-g e (o,o) 



1 



^ms 



0,0) 



(A.17) 



C- 



MS 



7 r-[G E (0,0)-G^(0,0)] 



1 



-/C£(0)-[G (0,0)-Gf s (0,0 



(A.18) 



C R (X) 

Here we used the fact that the difference [Ge (0,0) — G^ (0,0)] is independent of E. 
Rearranging eq. ( |A.18| ) using the fact that only the first two diagrams in the perturbative 
expansion of Gq(0, 0) are divergent yields 
1 1 



C Ws G! R (X) 



— lim 

r'-*0 



JC*(r>) + (r'\G° \r = 0) + (r'^KG^r = 0) 



+ 



(r' = 0\G° \r = 0W+ (r' = 0|G°KG°|r = 0W 



(A.19) 



Explicit calculation gives 

(r'|G°|r = 0) 



-M 



4nr' 

r'\G° V n G° \r = 0) = -°^^[\ - \n(m^r') - 7 + • • •] 



(A.20) 



where the ellipses represent terms that vanish as r' — > 0. In eq. ( |A.20| ) 7 is Euler's constant 
(7 ~ 0.577). Using dimensional regularization 



(r' = OlG^lr = 0)wm = 

U REG 



and 



(r' = 0|G8v;G8|r = 0)D 



IM 
REG 



-A-KarrM 2 I n (m 



(A.21) 



(A.22) 
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where I n (m n ) is the two loop integral 

T ( \ f ^ 



d n k 1 1 



(270™ J (27r)™ q 2 k 2 [(q-k) 2 + m 2 ] " 
Combining denominators with the Feynman trick 



f 1 f cf n q f d n k 1 
I n (m w ) = dx j (27r)nq2 J (^T [ k 2 + q 2 x(1 _ xj + TO 2 X ]2 



(A.23) 



dx 



d n q 7r n / 2 



T(2 - n/2) 



(A.24) 



(27r)™q 2 (27r) n [q 2 x(l - x) + m 2 x] 2 " n / 2 



Changing the momentum integration variable from q to p = yl — xq the above becomes 
r n/2 r c9_„ /9^ /-l ( _ _ _ _ r d n p 1 



(2tt 



p 2 [p 2 + m 2 ] 2 - n /2 



(A.25) 



Performing the p integration and the x integration gives 



7r n r(3 - n)(m 2 



2 Nn-3 



(27T) 



2/i 



(A.26) 



Consequently in the MS subtraction scheme 

(r' = 0|G°|r = 0)^ = 0, 
(r' = 0\G° %Gl\r = 0}- 



a T M 2 



MS 



8tt 



(A.27) 



Combining these results and using the small r behavior of /Cq(V) in eq. ( A. 3 ) yields 



1 



1 



Cms C r {\) 



and consequently, 



Gg*(0,0) = ^ 



8tt 



8tt 



I ^ ) 2,-1 



(^2) - ' 



(A.28) 



(A.29) 



The A dependence of 1/Cr(X) is exactly cancelled by the A dependence of In ( jjj appear- 
ing in the second term on the r.h.s of (|A.28| ) , leaving 1/Cjjg independent of A. Similarly, 
cancellation of the A dependent terms on the r.h.s. of ( |A.29|) leaves G^ s (0, 0) independent 
of A. Gg f5 (0, 0) is plotted in fig. 10 for a subtraction point of w = m n . 
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w 1 



50 100 150 200 250 300 

Pcm < MeV ) 



Fig. 10. The MS Green function at the origin normalized to f%, 

r E 



GM s (0, 0)/f%, as a function of |p| in MeV. The solid and dashed lines 



are the real and imaginary parts, respectively. 
Appendix B. Computing the effects of 2-derivative, 4-nucleon interactions 



In this appendix we show how to compute the sum of ladder diagrams in MS when 
a two-derivative, 4-nucleon operator is included. Consider the ladder sum including both 
the contact interactions ( 3.12 ) as well as the Yukawa part of one pion exchange: 



where V n is given in eq. ( |4.4[ ) as 
(P|K|P) = - 



v = v v + v e 



:b.v, 



(B.2) 



(p _p/ )2 + m 2> 8nf 2 

and f n = 132 MeV is the pion decay constant. (For the effective theory without pions, one 
can take the result we will derive and set = 0). The contact interaction V c is given by, 



(p\V c \p') = C 1 + 



P 2 + P /2 

2A 2 



(B.3) 



This can be conveniently rewritten as 

V c (p, p') =(7(1 + ME /A 2 ) - CM 



(E - p 2 /M) + (E — p' 2 /M) 
2A2 



(B.4) 



or in operator form as 

d 3 q d 3 q' 



V r = C 



(2tt) 3 (27r) : 



(1 + ME/A ) |q) (q'| - MiG^/A 2 , |q) (q' 



(B.5) 
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where (G^)" 1 = (E — Ho). 

The first term in eq. ( |rj.4| ) is easy to deal with, since E is a number, not an operator. 
Computing the effects of the second term involving the operator (G^) -1 requires some 
thought. There are three ways the (G'g) -1 term can enter the ladder diagram calculation. 

Firstly, there could be an insertion acting on the "external" legs (by external legs we 
mean nucleon propagators which interact via V n , but not through any contact interactions). 
This entails calculating the integral ((p| is an on-shell state) 



/ (£y^ (pi(^)" 1 ^(^)" 1 iq) = / (£yr (pl(^rHi + ^K)| q > 

= a f m, J (p| (G° E )~ 1 G E |q) , 



(B.6) 



where we made use of the fact that (p| {G Q E ) 1 |q) = for an on-shell state (p|; as well as 
the dimensionally regulated integral 

f ^ / 'if/ i \ f dn( l -4Tva n 
J (27r) n J (2yr) n (q - q') 2 + m 2 

This is equivalent making the replacement (G 1 ^) -1 — > a^m^. (Only the second factor of 
(G^) -1 in eq. (|B.6| ) comes from the interaction V c ; the first (Gg) -1 * s there to amputate 
the outgoing propagator). 

Secondly, one insertion of (G 1 ^) -1 cou ld act on internal nucleon lines (dressed by V n ). 
This gives rise to the integral 

II cPq dV ( , d {d0) -l ! /} = (q | (1 + GE V n) | q '> 

J J (2tt)" (2tt)" ^ hK E> lH/ J J (2tt)" (2tt)» V41 1 ^ ^ IH 7 

/y d n q d n q' - 
= o^m^ / / - — — - — — (q| Gk Iq') , 
J J (2n) n (27r) n x 1 1 7 

where we made use of the relation ( |2.8| ) between the full and free propagators, of the 
integral ( p.7| ) , and of the fact that J q 2r = in dimensional regularization. Again, 
(Gg) _1 just gets replaced by a^m^. 
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Finally, two insertions of (G E ) 1 could act on internal lines: 



= /•/• d^qd^q /q|(GO)-i(i + G g vr)|q/\ 
77 (27r)-(27r)- V4K EJ 1 B n) 7 

= a „^ / / U q U q ( q | (G%) G E |q'> 

r r d n q d n q' 
= K^) 2 yy (27r)n(2?r)n (q|G g |qO . 

Again we see that (G^) -1 § e ^ s replaced by a^m^. In conclusion, given eqs. ( p.4| ), (|B.5|) , 
the effect of including the 2-derivative operator in eq. ( |3.1| ) is simply to replace C by 

a^cf!-— + ™). (B,0) 

This is the result utilized in §4. The expression (|3.13| ) in §3 involves the substitution ( p.lOQ 
with a n set to zero (no pion contribution). 
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